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ABSTRACT 

We present an analysis of Spitzer-IRS observations of H2O, OH, HCN, C2H2 and CO2 emission, 
and Keck-NIRSPEC observations of CO emission, from a diverse sample of T Tauri and Hcrbig Ae/Be 
circumstellar disks. We find that detections and strengths of most mid-IR molecular emission features 
are correlated with each other, suggesting a common origin and similar excitation conditions for this 
mid-infrared line forest. Aside from the remarkable diff erences in molecular line st rengths between T 
Tauri, Herbig Ae/Be and transitional disks discussed in iPontoppidan et~a l. (2010fJ), we note that the 
line detection efficiency is anti-correlated with the 13/30 /im SED spectral slope, which is a measure 
of the degree of grain settling in the disk atmosphere. We also note a correlation between detection 
efficiency and Ha equivalent width, and tentatively with accretion rate, suggesting that accretional 
heating contributes to line excitation. If detected, H2O line fluxes are correlated with the mid-IR 
continuum flux, and other co-varying system parameters, such as L±. However, significant sample 
variation, especially in molecular line ratios, remains, and its origin has yet to be explained. LTE 
models of the H2O emission show that line strength is primarily related to the best-fit emitting area, 
and this accounts for most source-to-source variation in H2O emitted flux. Best- fit temperatures and 
column densities cover only a small range of parameter space, near ~ 10 18 cm -2 and 450 K for all 
sources, suggesting a high abundance of H2O in many planet- forming regions. Other molecules have a 
range of excitation temperatures from ~ 500— 1500 K, also consistent with an origin in planet-forming 
regions. We find molecular ratios relative to water of ~ 10~ 3 for all molecules, with the exception of 
CO, for which n(CO)/n(H20)^l. However, LTE fitting caveats and differences in the way thermo- 
chemical modeling results are reported make comparisons with such models difficult, and highlight 
the need for additional observations coupled with the use of line-generating radiative transfer codes. 



1. INTRODUCTION 

The chemistry of protoplanetary disks plays a pivotal 
role in the development of planetary systems, of hab- 
itable planets, and of life. However, our understand- 
ing of the chemical conditions in protoplanetary disks — 
especially in planet forming regions — is far from com- 
plete. In particular, there are many uncertainties about 
the role of volatiles such as water, which not only play 
a crucial role in life as we know it, but also in many 
other processes that affect the formation and evolution 
of planets. For example, the location of water-ice sub- 
limation, the so-called ice line, may define a boundary 
within which planetary cores do not grow large enough 
to become gas giants. And, molecular ices make up the 
bulk of many outer solar system bodies, including Kuiper 



Belt objects and outer-planet satellites. 

For many years, the detection of water in protoplane- 
tary disks re mained elusive, with only a few definitiv e ex- 
amples of ice (iMalfait et al.il998tlTerada et al 12007ft and 
vapo r (jNaiita et all 120001: ICaTr et al.l l2004flThi fc Bikl 
I2005D detections. Additionally, organic molecules 
were detected in a bsorption through a n edge-on disk 
Lahtiis et al.ll2006D. Re c ently, however, iCarr fc Naiital 
2008) and iSalvk et al.l (|2008l ) reported detections of 
a forest of emission lines of H2O, OH, HCN, C2H2 
and CO2 in high-resolution spectra from the Spitzer 
Space Telescope InfraRed Spectrograph (IRS), as well 
as of vibrational H2O lines with the Keck Near Infrared 
Spectrograph (NIRSPEC). Based on emitting temper- 
atures, fluxes and line shapes, these lines appeared to 
originate in the disk atmospheres at terrestrial-planet- 
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forming radii. Additionally, when present, the mid-IR 
water lines li kely dominate the cooling of the inner disk 
surface layer (|Pontoppidan et al"1l2010bl ). 

The frequent detections of water and other molecules 
in the inner regions of protoplanetary disks, made pos- 
sible largely by optimized high dynamic r ange observing: 
strat egies and data reduction routines (jCarr fc Naiital 

2008) , mark a turning point in our ability to study chem- 
ical evolution and molecular transport in the inner so- 
lar system. Because chemistry and line excitation are 
linked to disk physics, molecular emission can be used 
to study disk structure and physics. For example, wa- 
ter is destroyed by photodissociation and by high lev- 
els of ionization (e.g. iBergin et al.ll2003t iGlassgol d et al.l 

2009) , but is capable of self- shielding in certain condi- 
tions ([Bethell fc Bcrgin 2009]), and so water abundances 
may reflect the disk irradiation environment. Because 
water condenses at relatively small radii in the disk mid- 
plane (~ a few AU for solar- mass stars) , water vapor may 
be driven outward by vapor pressure differences, while 
water ice is carried inwards by planetesimals. Hence, 
water abundances may also reflect radial transport and 
planet formation timescales. Water abundances may 
similarly depend on vertical transport, which, when cou- 
pled with midplane condensation, coul d deplete water 
even well inside the canonical snow line (fMciicrin k et al.l 
2009). Finally, the strength of molecular emission lines 
depends upon the disk structure as a whole, including 
the gas temperature structure and the degree of dust 
settling/grain growth. 

One way to disentangle the many factors that affect 
molecular abundances and line strengths is to observe a 
large sample of disks, and look for trends predicted by 
the processes described above. Here wc report character- 
istics of a sample of disks observed with the Spitzer Space 
Telescope Infrare d Spectrograph (IRS), firs t described 
and presented in iPontoppidan et all (|2010bl ) (hereafter 
Paper I). In Paper I, we presented the overall character- 
istics of the spectra, including molecular detection rates, 
and reported a strong dependence of molecular emission 
on spectral type. In this paper, we analyze the IRS spec- 
tra in more detail, reporting correlations of line strengths 
and detections with system parameters, as well as molec- 
ular ratios calculated from LTE slab models. In addition, 
we use observations of CO rovibrational emission at 5 /im 
(G.A. Blake et al., in prep), to estimate molecular ratios 
with respect to CO. As one of the most abundant and 
easily-detected molecules in disks, CO is a tracer of the 
bulk of the disk gas, and provides a good baseline for 
comparing molecular abundances. M-band (^5 fim) CO 
rovibrational lines have similar excitation energies and 
collisional rates as the rotational H2O lines observed with 
the IRS, making these transitions particularly suited for 
comparison with the IRS observations. 

2. OBSERVATIONS AND REDUCTION 

2.1. Spitzer-IRS 

The mid-IR data in this study were all obtained with 
the Spitzer Infrared Spectrograph (IRS) in its high- 
resolution mode. The observed sample derives from a 
high S/N program designed for the detection of water 
and other molecular emission (PID 50641; PI J. Carr), 
as well as from a few programs in the Spitzer archive, and 



is described in detail in Paper I. A description of the ob- 
servations, including program IDs, integration times and 
Astronomical Observation Request (AOR) numbers, can 
also be found in Paper I. 

A brief explanation of the data reduction procedure is 
as follows. All data were reduced with our own IDL rou- 
tines , along with the IRSFRINGE package ( Tahuis et al.l 
2001. using a proce dure similar to that described in 
Carr fc Naiital (|2008[ ). We began with IRS droop- 
corrected, non-flat-fielded basic calibrated data (BCD) 
from the Spitzer Science Center IRS pipeline. In order 
to reduce noise caused by edge-effects, pipeline flat-fields 
were first divided by a low-order polynomial fit in the 
spectral and spatial directions. An off-frame was sub- 
tracted, data were flat-fielded, and noisy and bad pixels 
removed. Spectra were divided by an average of at least 
five standard stars, and minimally processed with the 
IRSFRINGE package. A more complete description of 
the reduction procedure can be found in Paper I. 

As reported in Paper I, the observed IRS spectra con- 
tain a forest of molecular emission lines from H2O, OH, 
HCN, C2H2 and CO2. For the analyses in Paper I and 
this work, positive detections require a 3.5<r peak (at 
two line locations for H2O and OH), with some line de- 
tections included or excluded through individual by-eye 
examination. Fluxes are calculated by defining and sub- 
tracting a linear fit to the local continuum, then sum- 
ming over the expected width of the line, except for the 
water emission lines, which are fit with Gaussians, as 
described in Paper I. Wavelength regions used are as fol- 
lows: OH (23.009-23.308 and 27.308-27.764 /im), HCN 
(13.837-14.075 (im), C 2 H 2 (13.553-13.764 /im) and C0 2 
(14.847-15.014 fim). All features overlap to some extent 
with H2O emission, with flux contributions from water as 
high as ~ 40% in the HCN and C2H2 regions, and even 
higher near the CO2 Q-branch. Thus, an H2O emission 
model (described in ST4J was subtracted before molecular 
line fluxes were computed. Nevertheless, this subtrac- 
tion does not have a significant effect on the statistical 
analyses we present. 

Statistical errors on individual points in the spectrum 
are calculated as a/\/N, where N is the number of frames 
used to calculate the flux, and a is the standard deviation 
of the values in that pixel. If N < 3, then the error is 
taken to be the standard deviation computed using three 
neighboring pixels in the dispersion direction. System- 
atic uncertainties are undoubtedly larger, especially for 
frames without dedicated background observations, but 
can be difficult to characterize. 

Line flux errors are taken to be the root squared sum 
of errors on individual points over the wavelength range 
used to compute the line flux. Upper limits are calculated 
assuming a 3.5cr line height (see Paper I for a discussion 
of this choice), and a feature shape defined by an LTE 
slab model, with column density, N, and temperature, 
T, taken to be near the mean of those derived in <gj 
However, upper limits are not very sensitive to the model 
parameters. Calculated fluxes and upper limits for all 
molecules except for water are shown in Tables Q] and [2j 
Detections and water fluxes can be found in Paper I. 

We also compute the flux in the 12.5-14 /jm interval 
(F13) and the 13/30 ^m SED slope (ni 3 _ 3 o), which are 
shown in Tables [4] and OH Precise definitions of these 
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parameters can be found in iFurlan et al.l ()2006[ ). Statis- 
tical errors are small; however, note that the variation 
in flux over the 12.5-14 /xm interval can be large (as 
much as several Jy), and that systematic errors can arise 
from background subtraction and the computation of the 
spectral response function. 

2.2. Keck-NIRSPEC 

M-band (~5 (im) spectra of CO Aw=l rovibrational 
emission were obtained for 54 of the 73 sources presented 
here. The data were obtained as part of a large M-band 
survey of Hcrbig Ac/Be (HAeBe) and classical T Tauri 
star (c TTs) disks, obtained from 2000-2010 with NIR- 
SPEC (jMcLean et al.l [20001 on the Keck II telescope. 
The spectra will be described in detail in G.A. Blake 
et al. (in preparation). However, we provide a few key 
results in this work, so as to compare with the Spitzer- 
IRS emission lines. Detailed ana lysis has been previously 
presented for transitional disks dSalvk et al.ll2008L I2009T ) 
and a few HAeBe disks (jBlake fc Boogertll2004[ ). repre- 
senting only a subset of the full survey. Line fluxes for 
the P(6) and P(12) transitions are shown in Tables [1] 
and [2j Upper limits are 3<r and assume a line width of 
5xl0" 4 urn (~ 30 km s" 1 ). 

3. CORRELATIONS 
3.1. Detection rate correlations between molecules 

Detections of most molecules are strongly correlated 
with each other, even when HAeBe and transitional disk 
non-detections (discussed in Paper I) are excluded from 
consideration. In Figure [T] we show the number of cTTs 
in each of four categories: both H2O and molecule X are 
detected, H 2 is detected but molecule X is not, molecule 
X is detected but H2O is not, and neither molecule is de- 
tected. For all molecules it is much more common for 
both or neither of the molecules to be detected. Addi- 
tionally, by-eye inspection of the data shows that many of 
the sources with only some molecular detections display 
evidence for additional molecules at lower S/N; therefore, 
it's likely that the various molecular emission features are 
even more strongly correlated. However, the C0 2 -only 
sources discussed in Paper I are a clear exception to this 
rule, and there may be others, such as Haro 1-4, which 
does not show strong evidence for H2O, but has strong 
HCN and C 2 H 2 features. 

In the lower left half of Table El we show probabilities 
(p) associated with linear regression of the detections, 
treated as binary variables (again, excluding HAeBe and 
transitional disks since they are all non-detections). For 
this and all future analyses, we define p in a standard 
way, as the probability of obtaining a linear regression 
slope at least as high as that observed, assuming the vari- 
ables are not related (i.e., the p- value associated with 
the t-statistic). Thus a probability p implies a corre- 
lation significant at the (1 — p) level. The small p for 
nearly all pairs of molecules shows that their detections 
are strongly correlated. A notable exception is C0 2 . 

3.2. Line flux and peak/continuum ratio correlations 
between molecules 

Line fluxes and peak/continuum ratios (a ratio of the 
sum of the continuum and maximum line flux density to 
the continuum flux density) are also strongly correlated 
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Fig. 1. — Number of cTTs for which H2O and molecule X are 
both detected, or only one is detected, or neither is detected. 
(HAeBe and transitional disks are excluded for clarity, since all 
are non-detections). 
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Fig. 2. — Peak/continuum correlations for OH (27.6 ^m feature), 
HCN, C2H2 and CO2 with H 2 (17.2 fim feature), along with best 
linear fits and associated probabilities. 

for several of the molecules. In Figure^ we show the OH, 
HCN, C 2 H 2 and C0 2 peak/continuum ratios against the 
H 2 17.2 /zm peak/continuum ratio, along with best-fit 
regression models and associated probabilities, and in 
Table [3] we show p for all molecular pairs. We see that 
peak/continuum ratios are correlated between several of 
the molecules, especially H 2 OH and HCN. Line fluxes 
(not normalized by the continuum) are correlated for all 
molecular pairs, but this is largely due to the fact that 
all line fluxes are correlated with continuum fluxes, as 
we will show. The fact that peak/continuum ratios are 
correlated for many molecules may mean that the line 
excitation is significantly influenced by qualities of the 
disk and radiation environment. 

Despite the correlations between peak/continuum ra- 
tios, variations in line ratios are still observed (and are 
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TABLE 1 
Line Fluxes: T Tauri Disks 
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3.66± 0.5 


3.37± 0.5 


V1121 Oph 


< 3.89 


< 3.52 


< 5.76 


< 3.79 


< 


3.02 


4.38± 0.5 


2.40± 0.4 


EC 82 


< 1.48 


< 1.44 


< 1.56 


< 1.06 


< 


0.88 


4.47± 0.1 


4.69± 0.1 



Note. — Fluxes arc 10~ 14 erg cm -2 s — . Upper limits are 3.5cr; error bars are 1 a. For CO lines, ns means the line is 
in emission but S/N or line coverage was not sufficient for determining a line flux, ea means the line is in emission, but a 
significant absorption component makes it difficult to determine the line flux, a means the line is in absorption. An ellipsis 
means the source was not observed at this wavelength. 



evident as the scatter in Figured]). In Figure El for exam- 
ple, we show two sources — DoAr 24E, in which the HCN 
and C2H2 have peak/continuum ratios at least twice that 
of neighboring H2O lines, and TW Cha, in which H2O 
features are comparable in strength to HCN and C2H2. 
Figure also shows a mean and scatter of the full sample 
(after normalization by the 17.2 /im line). Some features, 
including H2O and OH, appear to have little cross-sample 
variation, while HCN and C2H2 show significant varia- 
tion relative to H2O. 

3.3. Correlations between detections and system 
parameters 

In order to better understand why some sources show 
molecular emission and others do not, we searched for 



correlations between detections/non-detections and a 
number of stellar and disk parameters derived from the 
data or from the literature, including spectral type, X- 
ray luminosity, stellar mass, stellar luminosity, Ha equiv- 
alent width, accretion rate, disk mass, 13 /zm flux (F13) 
and SED slope (7^13-30 )■ The stellar and disk parame- 
ters used in this analysis are shown in Tables |4] and EJ 
and following Paper I, we consider the following disks to 
be transitional: CoKu Tau/4, HD 135344 B, SR 21, T 
Cha and LkHa 330. Spectral types and distances can be 
found in Paper I. 

In Figure [4j we show two major trends seen in the 
data — those between H 2 detections and (distance- 
normalized) F13 or ni3_3o. These trends are related 
to the conclusions presented in Paper I — that HAeBe 
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TABLE 2 

Line Fluxes: Herbig Ae/Be Disks 



Herbig Ae/Be disks 


Name 


OH 


OH 


HON 


C 2 H 2 


co 2 


CO 


(JO 




(23.2 /im) 


(27.6 pm) 








P(6) 


P(12) 


HD 36112 


<11.19 


<10.93 


< 8.43 


< 5.61 


< 4.52 


3.80± 0.2 


5.13± 0.1 


HD 244604 


< 2.99 


< 2.73 


< 4.68 


< 3.20 


< 2.44 


0.66± 0.1 


0.43± 0.1 


HD 36917 


< 3.71 


< 4.11 


< 5.72 


< 3.83 


< 2.82 


ns 


3.57± 0.6 


HD 37258 


< 2.36 


< 2.10 


< 4.97 


< 3.07 


< 2.50 


< 1.14 


< 0.23 


BF Ori 


< 2.28 


< 2.10 


< 3.50 


< 2.31 


< 1.75 


< 0.77 


< 0.20 


HD 37357 


< 3.41 


< 3.24 


< 4.23 


< 2.82 


< 2.20 


0.33± 0.0 


0.67± 0.1 


HD 37411 


< 3.03 


< 3.00 


< 2.02 


< 1.32 


< 1.05 


< 0.79 


< 0.28 


RR Tau 


< 2.75 


< 2.67 


< 3.86 


< 2.58 


< 1.95 


ns 


0.61± 0.2 


HD 37806 


< 8.64 


< 7.53 


<11.70 


< 7.98 


< 5.98 


1.68± 0.5 


2.80± 0.3 


HD 38087 


< 2.26 


< 2.53 


< 1.54 


< 1.08 


< 0.85 






HD 38120 


<13.87 


<12.79 


<11.37 


< 7.64 


< 6.22 


0.49± 0.0 


0.28± 0.0 


HD 50138 


<22.98 


<19.22 


<27.07 


<18.94 


<13.87 


6.72± 1.1 


13.86± 1.2 


HD 72106 


< 3.69 


< 3.48 


< 4.63 


< 3.14 


< 2.44 


< 0.85 


< 0.13 


HD 95881 


< 6.18 


< 5.52 


<10.25 


< 7.05 


< 5.22 






HD 98922 


<13.65 


<11.83 


<27.54 


<19.18 


<14.06 






HD 101412 


< 2.57 


< 2.39 


< 4.92 


< 3.26 


10.33 ± 0.55 






HD 144668 


<12.01 


<10.85 


<14.19 


< 9.67 


< 7.24 






HD 149914 


< 1.00 


< 1.05 


< 1.19 


< 0.81 


< 0.55 


< 0.81 


< 0.17 


HD 150193 


<13.84 


<12.21 


<13.20 


< 8.92 


< 6.93 


6.98± 0.4 


4.74± 0.3 


VV Ser 


< 3.24 


< 2.90 


< 5.37 


< 3.69 


< 2.77 


3.04± 0.3 


2.68± 0.3 


LkHa 348 


< 5.26 


< 4.49 


<11.24 


< 7.70 


< 5.60 


a 


a 


HD 163296 


< 9.01 


< 8.08 


<13.74 


< 9.38 


< 7.22 


8.50± 1.2 


9.06± 3.1 


HD 179218 


<19.55 


<18.59 


<16.96 


<11.68 


< 9.14 


ns 


2.33± 0.4 


HD 190073 


< 5.76 


< 5.14 


< 9.87 


< 6.56 


< 4.95 


4.15± 1.3 


3.00± 0.7 


LkHa 224 


< 5.37 


< 5.68 


< 8.66 


< 5.88 


< 4.50 


ns 


ns 



Note. — See Notes for Table Q] 



0.3 
0.2 

0.1 

0.0 
-0.1 
0.04 
0.03 
0.02 
0.01 

0.00 
-0.01 

0.10 
0.05 
0.00 



DoAr 24 E — fc — n - 














Normalized Population Mean and Scatter 




, II - -jv- r^J^' f\ 


L 3 njll-pj- ji - n Jl Jl 



13.0 



13.5 



14.0 14.5 

Wavelength [/U.m] 



15.0 



1 5.5 



16.0 



Fig. 3. — Top two panels show 13—16 /im spectra for two sources with significantly different molecular line flux ratios. Unlabeled 
features are primarily due to H 2 0. Below, we show mean (in black) and plus or minus one standard deviation (in gray) of the entire set 
of normalized spectra with observed H 2 emission. Hatched regions from left to right mark C 2 H 2 , HCN, OH, C0 2 and OH, respectively. 
Note the increase in population scatter at the locations of the HCN and C 2 H 2 lines. 



disks and transitional disks are much less likely to show 
emission. For example, there are no H2O detections for 
F13 > 6, which represent the HAeBe disks. Similarly, 
we see that there are no H2O detections for ni3_3o > 1, 
which corresponds to non-detections in transitional disks. 
However, Haro 1-16, which has detected H2O emis- 
sion, can be consi dered transitional by some definitions 
(jSalyk et al.l 120091 ) and shows evidence for inner clear- 
ing in millimeter maps ([Andrews et al.l[2009f ). making it 
a "borderline" transi tional, or perhaps pre-transitional 
(jEspaillat et al.ll2007t) . disk. 

To look for other interesting trends in the data with- 
out being influenced by the transitional and HAeBe non- 



detections, we also considered a restricted sample, which 
included only cTTs disks. Results are shown in Table [6l 
Several parameters show high statistical significance with 
detections; however, with 45 parameter/molecule pairs, 
we would expect a few pairs to be significant at the 95% 
level through chance alone, so these may or may not be 
significant. However, a notable result is that detections 
for all molecules are anti-correlated with ni3_3o at at 
least the 90% level. ^13-30 is a measure of the degree 
to which s mall grains have set tled out of the disk atmo- 
sphere (e.g. lFurlan et al.ll2006f) . with small ni3_3o corre- 
sponding to a higher degree of settling. Therefore, the 
observed correlation is consistent with greater detection 
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TABLE 3 

Classical T Tauri detection rate and peak/continuum correlation 

PROBABILITIES 



H 2 
OH 
HON 
C 2 H 2 
CO2 



H 2 

0.00 
0.00 
0.01 

0.12 



Peak / Continuum 
OH HCN C2H2 
0.00 0.00 
0.05 



0.00 
0.02 

0.14 



0.74 
0.20 
0.18 



0.00 

0.11 



co 2 

0.14 
0.20 
0.07 
0.04 



0.07 



NOTE — At lower left is the probability, p, associated 
with a linear regression of the two binary variables representing 
a detection or non-detection of the molecule. At upper right 
is the probability, p, associated with a linear regression of the 
peak/continuum ratios for each molecular pair. HAeBe and 
transitional disks were excluded. Bolded values have p < 0.05 and 
are statistically significant. 

rates in disks in which small grains have begun to settle 
to the midplane. This result follows naturally from the 
fact that settling of small grains allows us to see deeper 
into the disk, and hence observe a larger column and 
higher density (i.e., a density closer to the critical den- 
sity) of gas. However, the detailed (and likely complex) 
effects of grain settling on gas heating and structure have 
yet to be explored. 

We also find that all molecular detections are corre- 
lated with Ha at the 90% level or more, although, in- 
terestingly, only H2 O detections show a correlation with 
accretion rate. Water detections as a function of Ha and 
M are shown graphically in Figure [5] The correlation 
with M must be approached with some caution, as the 
M's derive from a variety of methods, which are known 
to have systematic offsets and significant scatter (e.g. 
iHerczeg &: Hill enbrand 2008]). In fact, since higher M's 
tend to be derived from UV excess measurements, while 
lower M's are derived from line observations, we find 
that reducing the UV-excess-derived values by factors of 
four to seven (the offset between UV-excess and Ha mea- 
surements found by Herczeg et al. 2008) can render the 
M detection trend statistically insignificant. However, 
the trend with Ha remains significant, and could mean 
that accretional heating is a necessary ingredient for the 
excitation of mid-IR molecular emission. 

3.4. Correlations between line fluxes and system 
parameters 

In Table [71 we show p for line fluxes versus various 
stellar and disk parameters. All molecular line fluxes 
are correlated with the continuum flux, F13, as well as 
with several other parameters. In Figure [51 we show the 
distance-normalized line flux of the 17.2 /zm H2O line 
against the 17 /im continuum flux. The dashed line has 
slope fixed at 1 , which is equivalent to assuming that the 
line flux varies in a sense proportional to continuum flux. 
We find that the H2O line flux is correlated with contin- 
uum flux, although the actual best- fit slope, 0.61 ± 0.14, 
differs somewhat from direct proportionality. Addition- 
ally, there is significant scatter about the trend, showing 
an intrinsic range in the line/continuum ratios among 
detected sources. Figure [5] also shows that HAeBe disks 
have line fluxes lower than the trend predicted by the 
cTTs emission lines, consistent with the conclusions of 
Paper I. 

As can be seen in Table [71 the correlation between 
continuum and line flux implies a correlation with several 
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Fig. 4. — H2O detections and non-detections as a function of 
distance-normalized continuum flux and ni3_3o color. Detected 
sources are shown by filled symbols. 
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Fig. 5. — H2O detections and non-detections in cTTs as a func- 
tion of Ha and M. Detected sources are shown by filled symbols. 

other system parameters that covary with the continuum 
flux. One such correlation, between line flux and L*, is 
shown in the second panel of Figure [6j and is marginally 
consistent with FocL. This could be related to increased 
heating of the gas, or an increase in the solid angle of 
the emitting region (due to gas being heated at larger 
distances from the star). In fact, such a relationship be- 
tween luminosity and solid angle has been found for CO 
(Pontoppidan et al. 2011; C. Salyk et al., in prepara- 
tion). The significance of the correlation observed here 
may, however, depend on a few stars with high mass ac- 
cretion rates. 

To look for other correlations with disk or stellar pa- 
rameters that are independent of the line-to-continuum 
correlation, we use peak/continuum ratios, thereby ac- 
counting for and removing the line/continuum correla- 
tion. However, we find that virtually all significant cor- 
relations with line fluxes are removed once we account for 
the continuum flux correlation. In fact, the number of 
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TABLE 4 
T Tauri System Characteristics 
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-0.451 ± 0.006 


Wa Oph 6 


35 




-6.64 




0.077 


39 


0.67 




I6I17I36I 


0.842 ib 0.002 


-0.386 ±0.005 


V1121 Oph 






-7.05 




0.028 


38 






m 


3.411 ±0.009 


-0.289 ±0.005 


EC 82 


5-11 
















[23] 


1.394 ±0.002 


0.475 ±0.003 



References. — (1) Alcala et al. 2000f (2) Andre & Montmerle 1994 (3) Andrews & Williams 2005 (4) Andrews & Williams 2007a (5) Andrews & 
Williams 2007b (6) Andrews & Williams 2009 (7) Bary, Weintraub & Kastner 2003 (8) Beckwith et al. 1990 (9) Brown ot al. 2007 (10) Brown et al. 

2008 (11) Casanova et al. 1995 (12) DAlessio et al. 2005 (13) Feigelson et al. 1993 (14) Fernandez et al. 1995 (15) Furlan ot al. 2006 (16) Garcia Lopez 
et al. 2006 (17) Gras- Velazquez, A., & Ray (18) Giidel & Telleschi 2007 (19) Hartmann et al. 1998 (20) Herczeg & Hillenbrand 2008 (21) Hughes et al. 
1994 (22) Isella, Carpenter & Sargent 2009 (23) Kessler-Silacci et al. 2006 (24) Kitamura et al. 2002 (25) Lommen et al. 2009 (26) Montesinos et al. 

2009 (27) Natta et al. 2006 (28) Neuhauscr et al. 1995 (29) Pinte et al. 2008 (30) Pontoppidan et al. 2008 (31) Pontoppidan et al., in prep (32) Prato, 
Greene & Simon 2003 (33) Robrade & Schmitt 2007 (34) Schegerer et al. 2006 (35) Seperuelo Duarte et al. 2008 (36) Shevchenko & Herbst 1998 (37) 
Spezzi et al. 2008 (38) Watson et al. 2009 

Note. — Error bars for F30 and 7113-30 are statistical errors. See the text for details. 

strong mid-IR molecular emission. (Again, it must be 
mentioned that Haro 1-16 has some characteristics of a 
transitional or pre-transitional disk, and has evidence for 
H2O emission). 

In Figure[8j we show that distance-normalized CO and 
H2O line fluxes are highly correlated (with p=0.004). In- 
terestingly, for most HAeBe disks, the observed H2O up- 
per limits are consistent with what would be predicted 
from the CO line fluxes. Thus, the data are at least 
consistent with no difference in CO/H2O line ratios be- 
tween cTTs and HAeBe disks, to within the constraints 



parameters with significant correlations is approximately 
what would be expected by chance alone. 

3.5. CO detections and correlations 

A majority of the disks in our sample were observed in 
the M band with Keck-NIRSPEC and most of these show 
CO rovibrational emission — including most HAeBe 
disks and many transitional disks. Figure [7] shows a 
histogram of detections and non-detections comparing 
CO and H2O. The majority of both HAeBe and tran- 
sitional disks show CO in emission, while neither show 
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Fig. 6. — Correlations between H2O line flux and distance-normalized continuum flux or stellar luminosity. Solid lines show best fits to 
sources with detected H2O emission. Dashed lines show best constant line/continuum fit and best fit with line flux proportional to stellar 
luminosity, respectively. Symbols are defined in Figure [4] 



TABLE 5 

HAeBe Continuum fluxes and colors 



Name 



' 13 



[Jy] 



"13-30 



HD 36112 
HD 244604 
HD 36917 
HD 37258 

BF Ori 
HD 37357 
HD 37411 
RR Tau 
HD 37806 
HD 38087 
HD 38120 
HD 50138 
HD 72106 
HD 95881 
HD 98922 
HD 101412 
HD 144668 
HD 149914 
HD 150193 

VV Ser 
LkHa 348 
HD 163296 
HD 179218 
HD 190073 
LkHa 224 



2.95 ± 0.01 
1.12 ± 0.01 
1.09 ±0.01 
0.95 ±0.01 

0.890 ± 0.006 
1.15 ±0.01 

0.663 ± 0.003 

1.242 ± 0.006 
7.30 ± 0.02 

0.064 ± 0.003 
6.22 ±0.02 
50.40 ±0.05 

I. 55 ± 0.01 
5.44 ± 0.02 

24.46 ± 0.05 
2.360 ± 0.008 

II. 35 ±0.02 
0.139 ±0.002 

9.50 ±0.02 
3.34 ±0.01 
6.82 ±0.02 

10.47 ± 0.02 
16.00 ±0.03 
4.57 ± 0.02 

8.0 ±0.2 



0.905 ±0.007 
-0.37 ±0.01 
0.34 ±0.01 
-0.91 ±0.02 
-0.97 ±0.02 
-0.16 ±0.01 
0.875 ± 0.008 
-0.20 ±0.01 
-1.034 ±0.006 

2.34±0.05 
0.243 ± 0.005 
-1.255 ± 0.003 

-0.26 ±0.01 
-1.391 ±0.008 
-1.310 ±0.004 
-1.007 ±0.007 
-1.037 ±0.006 

-0.50 ±0.03 
-0.534 ± 0.005 
-1.174 ±0.006 
-1.953 ±0.008 
-0.485 ± 0.004 
0.269 ± 0.003 
-1.34 ±0.01 
0.219 ±0.003 



Note. — See Table \4\ 

imposed by the upper limits on the H 2 line flux. In- 
stead, as shown in Figure HI we find that CO line-to- 
continuum ratios are systematically lower for HAeBe and 
transitional disks, just as was found for H 2 (Paper I), 
making the CO/H2O ratios similar for all classes of ob- 
jects. 

4. LTE MODELS 

4.1. Description of modeling procedures 

4.1.1. General modeling procedure 

To obtain a better physical understanding of the source 
of the disk emission, we fit the IRS data with emission 
models. In this work, we use an isothermal slab approx- 
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Fig. 7. — Histogram showing detections of mid-IR H2O and vi- 
brational CO emission. 



imation, in which the emitting area A, temperature T, 
and column density N are all free parame ters. These 
model s are the same as those utilized in iSalvk et al.l 
(2008), and assume a Gaussian local line shape with 
width (7=2 km s _1 . We make no correction for disk in- 
clination (since this has been measured for only a small 
subset of our sample), so the A we derive is a projected 
area, and N is a slant column density. 

The lines are not spectrally resolved with the Spitzer- 
IRS, and so the global line shape (which depends on 
large-scale, largely Keplerian motion, and disk orienta- 
tion) as well as the local line shape (which depends on 
the local gas velocity) are both unknown, but only the lo- 
cal line shape affects the slab model line flux. We choose 
a moderate value of the local line width, er, equivalent 
to the thermal velocity of H2 at ^700 K, or ~0.07i>Kep 
at 1 AU around a solar-mass star. This is a few times 
larger than the thermal velocity of H2O at 600 K and 
thus assumes some non-thermal gas motion. Since the 
actual local line broadening is essentially unknown, we 
consider here the effect of our chosen assumption. For 
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TABLE 6 

Regression probabilities for detections vs. cTTs parameters 



Molecule 


EW(Ha) 




log M 


M* 


M d 


i 




F13 


nl3.30 


H 2 


0.00 


0.54 


0.02 


(-) 0.21 


(-) 0.58 


0.21 


(-) 0.11 


0.16 


(-) 0.00 


OH 


0.04 


0.36 


0.23 


(-) 0.30 


(-) 0.79 


0.82 


(-) 0.20 


0.58 


(-) 0.00 


HCN 


0.02 


(-) 0.55 


0.28 


(-) 0.03 


0.86 


0.05 


(-) 0.05 


(-) 0.90 


(-) 0.06 


C2H2 


0.09 


(-) 0.86 


0.65 


(-) 0.10 


0.80 


0.04 


(-) 0.21 


(-) 0.97 


(-) 0.07 


CO2 


0.00 


0.30 


0.53 


(-) 0.50 


0.77 


0.85 


(-) 0.20 


0.42 


(-) 0.00 



Note. — Probabilities associated with the linear regression slope. Statistically significant correlations 
(p<0.05) are shown in bold. A (-) in front of the value means that it is an anti-correlation. 



TABLE 7 

Regression probabilities for line fluxes vs. cTTs parameters 



Molecule 


EW(Ha) 




M 


Mi 


M d 


i 


Li, 


F13 


nl3.30 


H 2 (15.2 nm) 


0.06 


0.04 


0.03 


0.00 


(-) 0.56 


-0.04 


0.00 


0.00 


(-) 0.63 


H 2 (17.2 fj.ro) 


0.08 


0.01 


0.02 


0.00 


(-) 0.46 


-0.05 


0.00 


0.00 


(-) 0.58 


OH (23.2 /im) 


0.55 


0.00 


0.05 


0.22 


0.81 


(-) 0.11 


0.05 


0.00 


(-) 0.60 


OH (27.6 /im) 


0.04 


0.12 


0.05 


0.02 


(-) 0.46 


(-) 0.07 


0.00 


0.00 


0.97 


HCN 


0.28 


0.92 


0.18 


0.21 


(-) 0.32 


(-) 0.21 


0.02 


0.03 


(-) 0.20 


C 2 H 2 


0.10 


(-) 0.83 


0.37 


0.29 


(-) 0.29 


(-) 0.12 


0.26 


0.00 


(-) 0.23 


CO2 


0.56 


0.03 


0.18 


0.01 


(-) 0.70 


(-) 0.19 


0.00 


0.00 


(-) 0.53 



Note. — Same notation as Tabic [B1 




-14 -13 -12 

Log P(6) Flux 

Fig. 8.— Correlations between CO P(6) and 17.2 ^tm H 2 line 
fluxes (in erg cm -2 s" 1 ) and M-band continuum flux density (in 
Jy), all normalized to 140 pc. Solid lines are best least-squares fits, 
while dashed lines are best constant line/continuum fits. Arrows 
denote 3<r upper limits. Filled symbols are sources with detected 
water, while open symbols are non-detections. 

molecules like H2O which have emission features at a va- 
riety of optical depths (r), the best-fit N scales roughly 
linearly with cr, since the data dictate r, and r oc N / a. 
For molecules like HCN, for which the emission is likely 
optically thin, the choice of a does not affect the total 
observed line flux, and thus does not affect the best-fit 
N, unless a is so small (< 0.1 km s _1 ), that the emission 
becomes optically thick. However, such small line widths 
are unlikely, as velocities cannot be sub-thermal, and T 
is several hundred K for all molecules. 

There are many limitations to using this simple slab 
model. The actual disk environment is quite complex, 
with densities and temperatures that depend on both 
height and radius, so the slab model only represents some 



average conditions in the disk atmosphere. Further, it 
has been shown that the line-excitation is not in LTE 
(jPontoppidan et all 12010a), as was predicted b y radia- 
tive transfer models (e.g. iMeiierink et al.ll2009l ). Many 
of the more complex issues related to line excitation and 
disk structure are being pursued with these sophisticated 
models. 

On the other hand, the slab models offer the benefit of 
providing general characteristics of the emitting gas for 
a large number of disks, with a relatively small number 
of complications and free parameters. Also, at least for 
H2O, non-LTE/L TE flux ratios remain c lose to unity for 
many transitions (jMeijerink et alJl2009f ). and the LTE 
results ap pear to reproduce the observed spectra reason- 
ably well (|Carr fc Naiitall2?M [Salvk et al.ll2008h . Here, 
we use the slab models to understand the bulk properties 
of the emitting atmospheres, to investigate the variety of 
emission observed, and to estimate molecular abundance 
ratios. 

4.1.2. Modeling procedure for H2 O 

The H2O emission spectrum observed by the IRS con- 
sists of lines with a large range of excitation temperatures 
and opacities, and so it is possible, in principle, to con- 
strain all three model parameters (A, T, and N). In par- 
ticular, A is principally determined from the strengths 
of the optically thick lines, while N and T both affect 
the relative strengths of lines as a fu n ction of wave- 
length. As discussed in iMeijerink et afl ()2009[ ). there is 
only a small change in line strength across a wide range 
of wavelengths, in contrast to predictions from simple 
disk models. Those authors interpret this observation as 
a depletion of water beyond ~ 1 AU, perhaps due to a 
"cold- finger" effect. The slab model cannot shed light on 
the radial distribution of the water abundance and the 
"flatness" of the spectrum is instead produced by a com- 
bination of high T and N. (However, the fact that the 
slab model provides as good a fit as it does is consistent 
with the "cold-finger" effect hypothesis, as the spectra 
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200 400 600 BOO 1000 1200 1400 1600 

Temperature [K] 

Fig. 9. — x 2 as a function of TV and T for VW Cha, along with 
10% and 5% probability contours. Each individual model is scaled 
by its own factor (A) to minimize residuals, and a is fixed such that 
Xr..,i = 1 f° r the best fit. White crosses show best- fit parameters 
used in our analysis, where for all molecules except for H2O, A is 
fixed to Ah 2 an d T is limited to lie between the vertical dashed 
lines (see also 34.1.31 . 

are reasonably well-fit by a single, relatively warm emit- 
ting temperature, equivalent to an inner disk emitting 
region with a sharp outer cutoff.) In Figure [3J we show 
the dependence of \ 2 on N and T, allowing A to be a 
tunable variable, and note that there is a well-defined 
minimum, but that N and T are somewhat degenerate. 

Because of the complexity of the water spectrum, and 
the limitations inherent to the simple model, no slab 
model provides a perfect fit to the data, and different fit- 
ting methods can produce different outcomes. We chose 
to compute 65 line peaks for data and model, and min- 
imize the difference between these two, with line peaks 
computed as the peak flux density minus an estimate of 
the local continuum. The set of lines includes any feature 



that could be reasonably distinguished from surrounding 
features, and the lines span wavelengths from 10-35 /im 
and thus a large range in excitation. This method does 
not require a global continuum fit, since it instead de- 
fines the continuum locally around each emission line. 
It is also not flux-weighted, and so takes weaker short- 
wavelength lines into account. We performed the fits by 
iterating over a grid of N and T, choosing a best-fit A for 
each of these pairs, and then choosing the best overall fit 
(based on the rms residuals). 

4.1.3. Modeling procedure for HCN, C2H2, CO2 and OH 

Continuum-subtracted spectra were used to fit all 
other molecular features. The continuum was deter- 
mined by fitting a spline to points along the spectrum 
where no strong H2O or other emission features were ex- 
pected. The preferred method and number of continuum 
points was carefully chosen to allow removal of broad- 
band features, such as those due to solid-state silicates, 
but not to remove narrower features due to gas molecules. 
However, uncertainties in continuum determination are 
unavoidable, and must be considered when determining 
the accuracy of any model fits. 

After continuum subtraction, the best-fit water model 
was also subtracted. We determined that our nominal 
H2O fits, which used the entire IRS wavelength range, 
did not produce a good match to the data below ~ 14/im. 
Therefore, we performed a separate fit to water lines in 
the 10-14 /jm region (which typically require tempera- 
tures increased by ~100 K and areas reduced by a factor 
of ~2), and subtracted these fits below 14 /mi. Finally, 
it was important to remove contributions from all other 
molecules; we began by ignoring these contributions, and 
then performed a few iterations of model fits and subtrac- 
tions. Figure [10] shows an example of the water model 
subtraction and its effect, for WX Cha. Some features 
are well-subtracted, but others are not, and these resid- 
uals affect the remaining molecular line shapes. 

For molecules except for water, the models are suffi- 
ciently degenerate that T, N and A can not be deter- 
mined simultaneously given our (primarily systematic) 
uncertainties. This is illustrated in Figure[9l which shows 
X 2 as a function of T and N, where A is adjusted to mini- 
mize residuals, and x 2 e d * s se * to 1 for the best fit. At low 
N (in the optically-thin regime) A and N are completely 
degenerate, and so \ 2 1S insensitive to changes in N . \ 2 
contours curve as N enters the optically-thick regime, as 
the line shapes begin to be affected by the opacity, but 
fits still remain somewhat degenerate with respect to all 
three fit parameters. Even when distinct minima exist 
(as in the case of CO2) they are not significant in a sta- 
tistical sense, because of uncertainties in the continuum- 
and water-model-subtracted data. Therefore, we have 
chosen not to fit all three parameters simultaneously. 

The following procedure was chosen to fit the molec- 
ular emission features. First, an appropriate range of T 
was determined by computing a normalized mean emis- 
sion feature and comparing with models, as shown in 
Figure [TT1 Then, fixing A to A^ 2 o, T (within the al- 
lowed range) and iV were varied to minimize % 2 (over 
a wavelength range wide enough to include all potential 
emission). A model with a different A could potentially 
produce a lower \ 2 (Figure [9]), but we believe that the 
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Fig. 10. — Example (for WX Cha) of H2O model removal and its effect on other molecular line shapes and strengths. The top panel 
shows the reduced data (in black) and best-fit H2O model (in gray). The bottom panel shows the data after H2O model removal, along 
with a model that includes HCN, C2H2 and C02- 
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Fig. 11. — Mean spectra, normalized to the line peak, in black. Gray areas delineate plus or minus 1 standard deviation from the mean. 
Normalized models, at temperatures meant to bracket the observed data, are shown in red and blue. 



systematic uncertainties do not justify the choice of this 
model over others. This diff ers from the procedure cho- 
sen by (|Carr fc Naiital 12008), who simultaneously fit all 
three model parameters for all molecules except C2H2. 
We stress that our modeling results are not inconsis- 
tent with theirs, but rather represent a more conservative 
approach based on our understanding of the systematic 
uncertainties. Some care must be taken with the inter- 
pretation of our results, as the effective emitting regions 
are almost certainly different for each molecule. How- 
ever, column densities can easily be adjusted as addi- 
tional data place constraints on A. The change in N 
would be directly anti-proportional to changes in A in 
the optically thin regime, but changes in N will be more 
significant if the lines turn out to be optically thick. 

To compute upper limits, we first estimated the emit- 
ting area, either by setting it equal to that derived for 
H2O (if observed) or by using the results of a linear fit 
to best-fit emitting area versus distance- normalized -F13. 
For each molecule, we then chose a characteristic tem- 
perature (shown as dashed lines in Figure IT2"j) and used 
the assumed emitting area to compute emission models. 
We used a 3c detection limit for the line (or band) peak 
to compute a column density upper limit. For H2O we 
used the 17.2 /xm line limits. 

4.1.4. Modeling procedure for CO 



We fit CO line fluxes for most CO-emitting sources 
with an LTE slab m o del fo llowing the procedure de- 
scribed in lSalvk et ail (|2009D . A few sources were omit- 
ted due to low S/N (HD 244604, HD 36917, RR Tau and 
HD 179218) and SR 21 was omitted beca use it is not 
well fit by an LTE model (|Salvk et al.ll2009r i. A rotation 
diagram was constructed from CO emission line fluxes, 
and these were fit with a grid of LTE models with T, N 
and A as free parameters. In addition, we used 13 CO 
fluxes or non-detections as a constraint on N. Finally, 
we restricted our analysis to those sources with at least 
two line flux measurements at J > 20, as we found the 
fits otherwise unreliable. 

4.2. Model Results 
4.2.1. H 2 Results 

Best-fit H2O model parameters and upper limits are 
shown in Tables [3] and O and a histogram of these re- 
sults is shown in Figure [T^l Best-fit T's cluster around 
450 K, and best-fit TV's cluster around 10 18 cm -2 . Tem- 
peratures of 450 K correspond to blackbody dust at a dis- 
tance of 1.65 AU at the surface of a passively heated disk 
with ?; = 4000 K and i?* = 2.5i? Q (jChiang fc Goldreichl 
1997). However, models predict that the gas in the 
disk's upper layers thermally decouples from the dust, 
so warm gas may exist at much larger distances from the 
star (e.g. lGlassgold fc Naiitdl2^0lTlKamp fc Dullemondl 
l2004f ). A^HsO is close to that derived from the water self- 
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Fig. 12. — Filled, gray histograms show best-fit model parame- 
ters from LTE fits. j4h 2 o is used for all molecules except for CO. 
Dotted, unfilled histograms show upper limits, while dashed verti- 
cal lines show the assumed temperature used to compute the upper 
limits. 

shielding process proposed by iBethell &: Berginl (|2009f) . 
but it should be kept in mind that the H 2 column den- 
sities in their models are derived in part from the IRS 
observations themselves. 

A key result of our fits is that the H 2 spectra are fit 
by a small range of N and T, despite the fact that line 
strengths vary greatly from source to source. The source- 
to-source variation in line strength is primarily explained 
by variations in A, with line strength and A being highly 
correlated. This is consistent with Figure [3j showing 
little variation in the H2O spectra after normalization. 

Our best-fit para meters are similar to those of 
iCarr fc Naiital (|2008f ). although some differences likely 
arise from their choice of wavelength region and fit min- 
imization. These choices have non-negligible effects on 
the model parameters, and there is no 'right' choice, since 
the model is much simpler than the true disk environ- 
ment. As an example, if we fit only the region shortward 
of 14 /im, we find slightly higher emitting temperatures 
(by ~ 100 K) and smaller emitting areas (by a factor of 
~ 2) than the nominal fits, consistent with the 12-14 /mi 
lines originating predominately in a disk region closer to 
the star than the longer wavelength lines. If we fit the 
continuum-subtracted spectrum, instead of the individu- 
ally measured line peaks, best-fit models differ from the 
nominal fit with higher T (by ~100 K) and N (by -0.5 
dex), and smaller A (b y ~4x). Our best -fit T differs 
more significantly from iSalvk et "all (|2008l) , because in 



that analysis A was fixed at a much smaller value, based 
on CO lineshapes; N remains similar for both analyses. 

4.2.2. HCN Results 

Figure fTTI shows the peak-normalized mean (after sub- 
traction of other molecular models) of all HCN-emitting 
spectra. The HCN lineshape is fit by a narrow range of 
temperatures — 600-800 K — because at low T the lo- 
cation of the line peak shifts longward, while at high T, 
the model develops a prominent hot-band "knee" on the 
short-wavelength side of the feature. 

The underlying H2O features contribute up to ~ 40% 
of the flux in this region. Therefore, in the case that 
the HCN is optically thin, the H2O model subtraction 
results in, at most, a 40% correction in N. However, if 
A << An 2 o then the emission may be optically thick, 
and uncertainties in the H2O flux subtraction could re- 
sult in larger uncertainties in N. With A = A^ 2 o, the 
best-fit model parameters are shown in Table [5] and Fig- 
ure [T2j Best-fit column densities range from a few xlO 14 
to a few xlO 15 cm~ 2 . 

There is an unaccounted-for residual in the ~13.8-13.9 
/itm region, apparent in Figure 111) but also visible in 
many individual spectra. One possible explanation is 
that the water has not been modeled correctly, but ob- 
served H 2 features do not line up well with the resid- 
ual. The residual does line up with the ,V2 = 5 — ^ 4 
and 1/2 = 4 — > 3 hot-bands of HCN, whose emitting lev- 
els are the lower states of the well-known 130 and 330 
/zm submillimeter masers (ICebbie et al.lUQfil . respec- 
tively. These high-lying bending states can be fed ei- 
ther by radiative excitation through combination bands 
involving the fairly weak C-N (v\) stretching mode, or 
though chemical reactions th at leave the produ ct HCN 
highly vibrationally excited (|Pichamuthul 119741 ) . Thus, 
the presence of this residual suggests that the bending 
vibrational levels of HCN may be non-thermally excited, 
and that caution must be used in the interpretation of 
the Q-branch profile using LTE models. 

4.2.3. C2H2 Results 

The mean C2H2 emission feature is shown in Figure fTTl 
after subtraction of both best-fit H2O models and best-fit 
HCN models. With A — Ah 2 o, the emission is optically 
thin, and the suite of emission features is consistent with 
T in the range 500-1200 K. Best-fit N range from 10 14 to 
10 15 cm~ 2 . The derived N will scale inversely with the 
assumed A until column densities of > 10 16 cm~ 2 when 
the emission starts to become optically thick. 

In addition to the residual at ~ 13.85 /im, three resid- 
ual peaks shortward of 13.6 /im are prominent in many 
sources, and their origin is unknown. At least one peak 
lines up with an H2O feature, and may be accounted 
for with a better H2O model. The 13.55 /im feature is 
at the location of a high-temperature OH emission line, 
and could be evi dence for highly-ex cited OH, like that 
seen m TW Hya (|Naiita et al.ll2010D . However, adjacent 
transitions in the OH rotational ladder are not appar- 
ent in the spectrum, and so this may not be the correct 
explanation. 

4.2.4. C0 2 Results 

The mean CO2 feature is shown in Figure [IT] The 
emission is mostly optically thin with A = Ah 2 q, though 
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TABLE 8 
Model Fit Results - cTTs Systems 



Name 


t h 2 o 


log Nh 2 o 


ToH 


log N H 


Thcn 


log Nhcn 


T C 2 H 2 


log Nc 2 H 2 


T co 2 


log N C o 2 


R 


LkHa 270 


450* 


<17.6 


1100* 


<15.3 


700* 


<15.4 


1200* 


<14.4 


700 


15.6 


1.30 


LkHa 271 


450* 


<17.4 


1100* 


<15.5 


700* 


<15.3 


1200* 


<14.0 


750* 


<15.1 


1.02 


LkHa 326 


650 


18.6 


900 


16.3 


600 


15.8 


1200* 


<15.0 


700 


16.3 


0.45 


LkHa 327 


350 


18.8 


1350 


14.9 


700* 


<15.2 


1200* 


<14.3 


800 


15.3 


2.4 


LkHa 330 


450* 


<17.4 


1100* 


<15.6 


700* 


<15.1 


1200* 


<14.2 


750* 


<14.9 


1.37 


LkCa 8 


450* 


<17.0 


1100* 


<14.9 


700* 


<14.9 


1200* 


<13.7 


750* 


<14.6 


0.990 


IQ Tau 


800 


17.6 


1100* 


<15.7 


800 


15.8 


1200 


15.6 


750* 


<15.5 


0.49 


V710 Tau 


450* 


<17.0 


1100* 


<14.8 


750 


15.0 


1200 


15.0 


800 


14.8 


1.00 


AA Tau 


400 


18.8 


1000 


15.6 


800 


15.6 


1100 


15.3 


600 


15.3 


0.74 


CoKu Tau/4 


450* 


<17.3 


1100* 


<14.9 


700* 


<15.0 


1200* 


<13.9 


750* 


<14.9 


0.980 


DN Tau 


450* 


<17.1 


1100* 


<14.9 


700* 


<14.9 


1200* 


<13.9 


750* 


<14.7 


1.04 


FX Tau 


450* 


<17.2 


1100* 


<14.9 


600 


15.0 


1200* 


<14.0 


750* 


<14.9 


1.02 


DR Tau 


500 


18.6 


1050 


15.6 


800 


15.6 


1200 


14.8 


400 


15.9 


1.1 


SX Cha 


350 


18.9 


1300 


15.3 


650 


14.7 


1200* 


<14.0 


750* 


<14.8 


1.0 


SY Cha 


450* 


<16.9 


1100* 


<14.6 


650 


15.0 


1200 


14.9 


800 


15.0 


1.04 


TW Cha 


450 


18.6 


1150 


15.3 


750 


15.4 


1200 


15.2 


750* 


<14.6 


0.76 


VW Cha 


500 


18.3 


900 


16.3 


800 


15.3 


1200 


14.8 


400 


15.8 


1.2 


VZ Cha 


450 


18.6 


1100 


15.8 


800 


15.5 


800 


15.4 


750* 


<14.9 


0.78 


WX Cha 


450 


18.8 


950 


15.6 


800 


15.6 


1200 


14.7 


400 


15.7 


0.79 


XX Cha 


250 


20.3 


700 


16.0 


800 


15.0 


1200 


15.7 


700 


15.2 


0.88 


T Cha 


450* 


<16.6 


1100* 


<14.6 


700* 


<14.4 


1200* 


<13.4 


750* 


<14.2 


0.960 


Sz 50 


450* 


<17.0 


1100* 


<14.7 


700* 


<14.7 


1200* 


<13.9 


750 


15.0 


1.04 


HD 135344 B 


450* 


<16.8 


1100* 


<14.8 


700* 


<14.6 


1200* 


<13.6 


750* 


<14.3 


1.03 


HT Lup 


450* 


<17.7 


1100* 


<15.0 


700* 


<15.4 


1200* 


<14.5 


800 


15.4 


1.57 


GW Lup 


450* 


<17.1 


1100* 


<14.9 


700* 


<14.9 


1200* 


<13.8 


750 


15.2 


1.01 


GQ Lup 


400 


17.9 


950 


15.3 


600 


15.2 


1200* 


<13.8 


750* 


<14.3 


1.2 


IM Lup 


450* 


<17.1 


1100* 


<14.8 


700* 


<14.8 


1200* 


<14.1 


750 


15.0 


1.11 


HD 142527 


450* 


<17.4 


1100* 


<15.5 


700* 


<15.1 


1200* 


<14.3 


750* 


<14.9 


2.97 


RU Lup 


500 


19.0 


1100* 


<15.2 


800 


15.7 


1200* 


<14.5 


400 


16.3 


0.90 


RY Lup 


450* 


<17.5 


1100* 


<15.2 


700* 


<15.2 


1200* 


<14.2 


750* 


<15.0 


1.21 


EX Lup 


300 


19.0 


1600 


15.3 


700* 


<14.5 


1200* 


<13.7 


750* 


<14.3 


2.0 


AS 205 


300 


20.6 


700 


16.8 


800 


15.2 


1200 


15.3 


600 


15.6 


2.1 


Haro 1-1 


450* 


<17.0 


1100* 


<15.0 


700* 


<14.9 


1200* 


<13.7 


750* 


<14.6 


0.980 


Haro 1-4 


450* 


<17.2 


1100* 


<15.0 


700 


15.0 


1200 


15.0 


750* 


<14.8 


1.05 


VSSG1 


350 


17.8 


1100* 


<14.6 


600 


15.0 


950 


15.6 


800 


14.6 


2.1 


DoAr 24E 


450 


18.3 


1100 


15.3 


800 


15.8 


1050 


15.9 


700 


15.4 


1.2 


DoAr 25 


450* 


<16.9 


1100* 


<14.9 


800 


15.1 


1200 


14.8 


750* 


<14.6 


1.00 


SR 21 


450* 


<17.4 


1100* 


<15.4 


700* 


<15.1 


1200* 


<14.1 


750* 


<14.9 


1.43 


SR 9 


450* 


<17.4 


1100* 


<14.9 


700* 


<15.1 


1200* 


<14.1 


750* 


<14.9 


1.07 


V853 Oph 


300 


18.6 


1300 


14.6 


700 


14.7 


1050 


16.0 


750* 


<14.5 


1.2 


ROX 42C 


450* 


<17.2 


1100* 


<15.0 


700* 


<15.0 


1200* 


<13.9 


750* 


<14.8 


1.05 


ROX 43 A 


450* 


<17.4 


1100* 


<14.8 


700* 


<15.1 


1200* 


<14.1 


750* 


<14.9 


1.36 


Haro 1-16 


450 


18.3 


1200 


15.3 


800 


15.1 


1200* 


<14.1 


750* 


<14.7 


0.79 


Haro 1-17 


450* 


<16.8 


1100* 


<14.8 


700* 


<14.7 


1200* 


<13.4 


750* 


<14.5 


0.980 


RNO 90 


450 


18.3 


1000 


15.6 


750 


15.5 


1200 


14.7 


800 


15.2 


1.8 


Wa Oph 6 


400 


19.3 


800 


16.6 


700 


15.3 


1200* 


<14.4 


500 


16.0 


0.60 


V1121 Oph 


250 


20.9 


1100* 


<14.8 


700* 


<14.9 


1200* 


<14.0 


750* 


<14.7 


1.7 


EC 82 


450* 


<17.2 


1100* 


<15.0 


700* 


<14.9 


1200* 


<14.1 


750* 


<14.7 


1.95 



NOTE. — Units are cm 2 for N and AU for R (where R = y/A/w). 
Denotes a fixed temperature for computing upper limits on N . 



slightly optically thick at the center of the Q-branch 
peak. The complete suite of feature shapes is consistent 
with T from 500-1200 K, but best-fit T cluster near 700 
K. Best-fit N are ~ 10 14 — 10 16 cm -2 and are inversely 
proportional to A until N > 10 17 cm~ 2 . 

No significant residuals are noted; those seen in Figure 
[TT1 are most likely H2O model subtraction residuals or 
data-reduction artifacts. 

4.2.5. OH Results 

OH was fit with emission doublets from 17 to 30.5 ^m. 
Line strengths decrease towards shorter wavelengths. 



Derived excitation temperatures are strongly dependent 
on the strength of the relatively weak short wavelength 
lines, which are in turn dependent on the continuum sub- 
traction, making T somewhat more uncertain for OH 
than for other molecules. In addition, with A = An 2 o, 
the low energy lines are beginning to become optically 
thick, and so the shape of the emission spectrum is in- 
fluenced by both N and T. 

Figure [T3l shows the rather large range of T (700-1600 
K) that may be appropriate for fitting OH, with the mean 
profile in this case composed of spectra normalized to 
the 30.7 fj,m line peak. A much higher T is inconsis- 
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TABLE 9 

Model fit results - HAeBe Systems 



Name 


T H 2 


log Nh 2 o 


ToH 


log N H 


Thcn 


log Nhcn 


T C 2 H 2 


log Nc 2 h 2 


T co 2 


log N C o 2 


R 


HD 36112 


450* 


<17.7 


1100* 


<15.6 


700* 


<15.3 


1200* 


<14.4 


750* 


<15.1 


2.13 


HD 244604 


450* 


<17.8 


1100* 


<15.4 


700* 


<15.5 


1200* 


<14.6 


750* 


<15.3 


2.54 


HD 36917 


450* 


<17.9 


1100* 


<15.7 


700* 


<15.6 


1200* 


<14.7 


750* 


<15.4 


2.51 


HD 37258 


450* 


<17.9 


1100* 


<15.4 


700* 


<15.6 


1200* 


<14.9 


750* 


<15.4 


2.90 


BF Ori 


450* 


<17.7 


1100* 


<15.4 


700* 


<15.5 


1200* 


<14.6 


750* 


<15.3 


2.30 


HD 37357 


450* 


<17.8 


1100* 


<15.6 


700* 


<15.5 


1200* 


<14.6 


750* 


<15.3 


3.17 


HD 37411 


450* 


<17.7 


1100* 


<15.8 


700* 


<15.4 


1200* 


<14.4 


750* 


<15.2 


2.48 


RR Tau 


450* 


<17.7 


1100* 


<15.4 


700* 


<15.5 


1200* 


<14.5 


750* 


<15.2 


13.1 


HD 37806 


450* 


<17.4 


1100* 


<15.1 


700* 


<15.2 


1200* 


<14.3 


750* 


<14.9 


7.19 


HD 38087 


450* 


<18.2 


1100* 


<16.0 


700* 


<15.9 


1200* 


<14.8 


750* 


<15.8 


1.15 


HD 38120 


450* 


<17.6 


1100* 


<15.4 


700* 


<15.2 


1200* 


<14.4 


750* 


<15.0 


7.10 


HD 50138 


450* 


<16.9 


1100* 


<14.6 


700* 


<14.7 


1200* 


<13.9 


750* 


<14.4 


11.5 


HD 72106 


450* 


<17.7 


1100* 


<15.4 


700* 


<15.4 


1200* 


<14.4 


750* 


<15.1 


2.22 


HD 95881 


450* 


<17.3 


1100* 


<14.9 


700* 


<15.1 


1200* 


<14.2 


750* 


<14.9 


1.80 


HD 98922 


450* 


<17.2 


1100* 


<14.8 


700* 


<15.0 


1200* 


<14.2 


750* 


<14.8 


27.5 


TTF1 1 01 Al 9 






1 1 on* 




700* 


^ 1 
SlO.l 


1 900* 




soo 


1 « 
1D.U 


l.DD 


HD 144668 


450* 


<17.3 


1100* 


<15.0 


700* 


<15.0 


1200* 


<14.2 


750* 


<14.8 


4.02 


HD 149914 


450* 


<17.1 


1100* 


<15.0 


700* 


<14.8 


1200* 


<14.3 


750* 


<14.7 


1.00 


HD 150193 


450* 


<17.4 


1100* 


<15.1 


700* 


<15.1 


1200* 


<14.2 


750* 


<14.8 


2.74 


VV Ser 


450* 


<17.4 


1100* 


<15.0 


700* 


<15.1 


1200* 


<14.3 


750* 


<14.9 


4.33 


LkHa 348 


450* 


<17.4 


1100* 


<14.9 


700* 


<15.2 


1200* 


<14.3 


750* 


<14.9 


6.12 


HD 163296 


450* 


<17.3 


1100* 


<14.9 


700* 


<15.0 


1200* 


<14.1 


750* 


<14.8 


2.39 


HD 179218 


450* 


<17.3 


1100* 


<15.1 


700* 


<15.0 


1200* 


<14.1 


750* 


<14.8 


5.43 


HD 190073 


450* 


<17.5 


1100* 


<15.1 


700* 


<15.3 


1200* 


<14.4 


750* 


<15.1 


9.19 


LkHa 224 


450* 


<17.2 


1100* 


<14.9 


700* 


<15.0 


1200* 


<14.1 


750* 


<14.8 


15.5 


Note. — Units are 


cm -2 for N ar 


id AU for R (where R 


= VAJ 


iF). 













Denotes a fixed temperature for computing upper limits on N . 

tent with the data, suggesting a difference between the 
emission observed here primarily in the Long-High mod- 
ule, and the OH emission observed in th e Short-High 
modu le in the transitional disk TW Hya ([Naiita et al.l 
120101) . though the latter preferentially probe lines with 
higher excitation energy. Potential evidence for higher 
excitation lines in our data is the residual line at 13.55 
(im (discussed in £14.2.31) which is coincident with a high- 
excitation OH transition. High excitation temperatures 
may be indicative of non-thermal excitation, such as in 
the c ase of the "promp t" emission observed in comets 
(e.g. iBonev et al.l 1200a ). in which rotationally-excited 
OH is produced via the photolysis of H 2 0. For the most 
part, the moderate excitation temperatures we observe 
are instead reasonable thermal temperatures for few-AU 
disk atmospheres around cTTs disks. 

Best-fit temperatures are >800 K and typically lie be- 
tween 900 and 1100 K. Best-fit N are near 10 fs cm" 2 . 
The strongest lines are already becoming optically thick 
at these column densities, so N could differ significantly 
if A 7^ Ah 2 o- For example, in a test with the emitting 
area reduced by a factor of 10, we find the der ived N to be 
increa sed by about a factor of 25. Similarly. iSalvk et al.l 

(2008) found column densities of 2 x 10 1 ' cm -2 when 
fitting 3 fim OH emission features for DR Tau and AS 
205, which are a factor of ~50 higher than those found in 
this analysis; however, emitting areas were fixed to values 
27 and 30 times smaller, respectively. These N are sig- 
nifican tly smaller than those found by Beth ell fc Berginl 

(2009) to be consistent with water self-shielding, but we 
note that higher N of the right order for this process are 
not ruled out with smaller A. 



4.2.6. CO Results 

Best-fit model parameters can be seen in Table fTOl and 
Figure IT21 Temperatures are typically > 900 K, with a 
cluster at very high T (> 1500 K), though temperatures 
this high may be unrealistic. Despite the large range 
of T, N ranges over only about 1 dex, from ~10 18 — 
10 19 cm -2 . 

4.3. Molecular Ratios relative to -ff 2 O 

Figure [14] shows molecular ratios with respect to 
H 2 0. N x /Nn 2 o~ 10~ 3 for all molecules. Also 
plotted are ratios derived fro m previous disk ob- 
servations (ICarr fc Naiital [20 081 , ice observations in 
quiesce nt clouds dWhittet etall [20071) . high -mass 
YSOs dGerakines et al.l [1999D and low-mass YSOs 
(jPontoppidan et al.l 120081), as well as severa l disk 
models dMarkwick et all l2002t lAgundez et al.1 120081: 
Willa cy fc Woods! I2009D . Additionally. iGlassgold e4~atl 
(2009) (not plotted) predict OH/H 2 ratios in the warm 
molecular layer of the disk similar to the ratios found 
here — 10~ 4 — 10~ 2 for a variety of chemical models, 
sampled at radii from 0.25-3 AU. 

The large diffe r ences between our work and that 
of ICarr fe Naiital ([2001 arise from differences in de- 
rived column densities for both H 2 and any partic- 
ular molecule, which are due to differences in choice 
of fitting method and assumptions, as discussed above. 
With Spitzer-IRS, the Q-branch transitions are com- 
pletely blended together, but at higher spectral resolu- 
tion individual Q, P, or R-branch lines could be resolved, 
which would make possible independent determinations 
of emitting area. This highlights the need for higher- 
spectral-resolution and, eventually spatial resolution, ob- 
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Fig. 13. — Mean (black solid line) plus or minus the standard deviation (gray) of normalized OH emission sources, after subtraction of 
other molecular emission models. Plotted in red, blue, and dotted green, respectively, are optically thin 700, 1600, and 4000 K models, 
normalized to the 30.5 fim peaks. 



TABLE 10 
CO Model Fit Results 



Name 


Tco 


log N C o 


Rco 


LkHa 326 


1550 


19.3 


0.08 


LkHa 327 


1625 


19.5 


0.10 


LkHa 330 


850 


17.4 


0.30 


LkCa 8 


1475 


14.6 


2.41 


V710 Tau 


1675 


19.3 


0.05 


AA Tau 


950 


18.6 


0.13 


DR Tau 


1250 


18.9 


0.17 


HD 135344 B 


900 


17.9 


0.10 


GQ Lup 


1675 


18.6 


0.08 


AS 205 


1000 


18.7 


0.35 


Haro 1-4 


1675 


18.5 


0.04 


VSSG1 


1600 


18.6 


0.09 


DoAr 24E 


1275 


18.3 


0.21 


SR 9 


1575 


16.9 


0.16 


V853 Oph 


975 


18.0 


0.09 


Haro 1-16 


1250 


18.2 


0.08 


RNO 90 


1625 


18.2 


0.22 


Wa Oph 6 


1675 


18.5 


0.10 


V1121 Oph 


1125 


18.5 


0.12 


EC 82 


1000 


17.9 


0.39 


HD 36112 


775 


17.9 


0.51 


HD 37357 


575 


19.5 


0.61 


HD 37806 


625 


19.8 


1.39 


HD 38120 


800 


18.3 


0.27 


HD 50138 


725 


19.1 


0.97 


HD 150193 


725 


18.5 


0.47 


VV Ser 


1175 


18.1 


0.61 


HD 163296 


850 


18.5 


0.32 


HD 190073 


525 


20.4 


2.64 


Note. — Units are cm 


~ 2 for N 


and AU 



for R (whore R = yjApK). 

servations of molecular emission. The brightest sources 
are being pursu ed at high resolution with c urrent instru- 
mentation (e.g. IPontoppida n et aL1l2010al ). but a more 
complete census of water-emitting disks may require the 
next generation of ground- and space-based mid-IR spec- 
trographs. 

C0 2 /H 2 ratios in ices have been observed to be 
~ 0.1 — 1 in a variety of en v ironments dGerakines et al.l 
Il999t iWhittet et alj 120071: IPontoppidan et al.l l2008fh 
higher than our obse rved values and at th e upper range 
of those derived by ICarr fc Naiital (|2008[ ) (see the far 
right of Figure lUt)) . This may reflect chemical evolution 



towards lower values in th e inner disk region, as s een, for 
example, in the models of iMarkwick et all (|2002l ). 

Figure [TJ] also shows molecular ratios derived from 
thermo-chemical disk models. Disk-integrated c olumn 
density ratios at 1 AU from Markwi ck et al.l (|2002l ) 



(circles) and ([Willacy fc Woods! [2009h (stars) are plot- 
ted to the right side of the observations, while column 
density ratios from the disk surface (T>150 K; stars; 
iWillacy h Woods! 120091 ) a nd from photo-domin ated re- 
gions (PDR) modeled by lAgundez et~aTl (|2008| ) (trian- 
gles) are plotted to the left. Although we show ratios 
based on integrated vertical columns, these are likely not 
relevant for the disk surface being probed by the IRS 
molecular emission. 

It is immediately apparent that comparisons of data 
and models can be quite complex. For example, orders of 
magnitude differences in molecular ratios are computed 
when considering different regions of the disk, both radi- 
ally and vertically, for many molecules. And, the emitted 
lines themselves originate from a range of disk radii and 
heights, which are not ev en the same from tr ansition to 
transition. For example, iWoitke et al.l (|2009f ) find that 
pure rotational H2O lines could in principle originate in 
both the inner and outer (>10 AU) disk, and even in- 
side the disk inner rim. Nevertheless, the scatter ob- 
served in column density ratios is smaller than some of 
the model differences, holding promise for observations 
to eventually distinguish between models. One thing we 
can already say is that the observed abundances of C2H2 
re quire active car b on ch emistry, such as that described 
in lAgundez et al.l (|2008l ). as models without such net- 
works produce nowhere near e nough C2H2 in th e inner 
disk atmosphere; for example, Will acv fc Woods! (2009) 
find C 2 H 2 /H 2 ratios of ~ 10" 10 (not plotted). Alterna- 
tively, the high abundance of C2H 2 could be produced by 
the burning of polycyclic aromatic hydrocarbons (PAHs; 
though the production of PAHs may require C 2 H2 as a 
precursor; see discussion in ^6.4[) . 

A more promising approach to comparing models and 
observations may be to couple thermo-chemical mod- 
els with line-generating radiative transfer code, and 
directly compare output model spectra with observa- 
tions. This approa ch has been pursued, for example, by 
iKamp et al.l ((2010) for atomic fine-str ucture lines and 
low-J rotati onal CO line s, as w ell as [Meiierin k et al.l 
(|200l) and IWoitke et all (|20ol ) for far-IR rotational 
lines of H2O. An alternative approach is to map abun- 
dance structures using observations and radiative trans- 
fer, and compare these to the abundances predicted by 
various thermo-chemical models. This approach has al- 
ready yielded interesti ng results about water abundances 
(|Meiierink et al.|[2009D , and we expect that it will lead to 
significant progress in the comparison of chemical models 
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Fig. 14. — Molecular ratios, relative to H2O (filled squares), with 
error bars showing the Icr scatter . Also shown are p revious ob- 
servations in disk s (open square s: Carr & Naiita 2008), quiescent 
cloud ices (cross: Whittet ct al. 2007), high-mass YSO ices (open 
diamond: Gcrakincs ct al. 1999) a nd low -mass YSO ices (open in- 
verted triangle: [Pontoppidan ct al. 20 03) ■ To the left are models of 
the disk surface (open stars: Willacv & Woods 2009) and photon- 
dominated regions (open triangles: Aeiindcz ct al. 2008); to the 
right are vertically integrated model values at 1 AU (open stars: 
Willacy & Woods 2009; circles: Markwick et al. 2002). 

and observations in the near future. 

4.4. Molecular ratios relative to CO 

In Figure flS) we show molecular ratios with respect to 
CO. We find H 2 0/CO rat ios slightly above 1 — nearly 
the same as that found by lCarr fc Najital ()2008l ) and ly- 
ing bet ween vertically-i n tegra ted and disk surface model 
ratios. iGlassgold et al.l (|2009l ) also find H2O/CO ratios 
near 1 for several types of models. Again, other molecu- 
lar rat ios are different from those found in lCarr feTN aiita 
(2008) because of our choice of modeling assumptions. 
The scatter in H2O/CO ratio is only ^0.5-1 dex and so 
these ratios may prove helpful in distinguishing between 
chemical models. Note, however, that we have found sig- 
nificantly different emitting areas for CO and H2O and 
so the NIRSPEC and IRS emission lines are almost cer- 
tainly probing different emitting regions. 

5. SUMMARY OF RESULTS 

We have presented and analyzed Spitzer-IRS spectra 
of H 2 0, OH, HCN, C 2 H 2 and C0 2 and Keck-NIRSPEC 
spectra of CO from a large sample of cTT and HAeBe 
disks. In addition to the detection dependence on spec- 
tral type and SED class (transitional versus classical) 
noted in Paper I, we also find that detection efficiency 
depends on the disk color, rti3_3o, Ha equivalent width 
and, tentatively, on M. Radiative transfer disk mod- 
els suggest a number of ways to enhance line emission, 
including increasing the molecular abundance, lowering 
the dust /gas rat io in the disk atmosp here, and increas- 
ing gas heating (jMeijerink et al.l 120091 ). Since 7113-30 is 
a measure of the settling of small grains in the disk at- 
mosphere, these data suggest that it may be necessary 
to lower the dust/gas ratio in the atmosphere in order 
to produce these lines. In addition, the correlation with 




H 2 OH HCN C 2 H 2 C0 2 
X 



Fig. 15. — Molecular ratios, rela tive to CO. References and sym- 
bols are the same as in Figure fl4l 

Ha suggests that accretion is also important, perhaps 
because it provides additional heating for the excitation 
of lines. 

We find that detections, line fluxes and line-to- 
continuum ratios for nearly all molecules are correlated, 
suggesting a common origin or common excitation condi- 
tions for this forest of lines. Line fluxes (when detected) 
are also correlated with continuum flux and many other 
covarying system parameters (including L*, M* and M), 
and once this variation is accounted for, no other statisti- 
cally significant correlations are noted between line fluxes 
and system parameters. The observed correlation could 
be the result of any of these physical parameters, but a 
good candidate is the stellar luminosity, which provides 
heating for the line emission, and will also determine the 
size of the line-emitting region. 

Although line fluxes are correlated, line flux ratios (for 
example, HCN/H2O) do vary from source to source. No 
obvious trends in line flux ratios are seen, but the varia- 
tions suggest real chemical differences whose origin may 
be revealed with future study. 

We also compare the mid-IR emission lines with 
near-IR CO vibrational emission observed with Keck- 
NIRSPEC. Interestingly, CO lines are often seen even 
when mid-IR emission is not, including in many HAeBe 
and transitional disks. However, rather than find evi- 
dence for different H 2 0/CO ratios, we instead find that 
both transitional disks and HAeBe disks have system- 
atically lower CO line/continuum ratios, consistent with 
the reduced IRS line/continuum ratios reported in Paper 
I. This suggests possible systematic differences in line ex- 
citation between the two types of sources. 

While H2O line strengths do differ significantly from 
source to source, the emission spectra are similar if nor- 
malized by line strength. In the context of an LTE model, 
the line strength is primarily determined by the emit- 
ting area, and so this accounts for most of the variation. 
Temperature and column density are then determined by 
the spectral slope. Best-fit temperatures lie near 450 K 
and N~ 10 18 cm -2 for all sources, though results depend 
somewhat on the choice of fitting method and wavelength 
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region. The high temperatures suggest an origin in the 
planet-forming (few AU) region of the disk, and TV is 
consistent with chem ical models of H2O formation (e.g. 
IGlassgolde t al. 200H) as well as with m odels that include 
UV d estruction and H2O self-shielding (jBethell fc Berginl 

Other molecules are fit with a range of excitation 
temperatures, but all are representative of the planet- 
forming region of the disk, some 500 K or higher. Be- 
cause of model degeneracies, coupled with systematic un- 
certainties, we chose to fix emitting areas, but emitting 
areas and column densities can in the future be scaled 
to match results from additional observations. With this 
technique, we derive molecular ratios with respect to wa- 
ter and CO of - 10~ 3 for HCN, C 2 H 2 , C0 2 and OH. 
We find CO/H2O rati os near 1, similar to values from 
previous obs ervations (|Carr fe Naiitall2008[) and several 
disk models (iMarkwick et al.ll2002b lAeundez et al.ll2508t 
iGlasseold et alj|200ft iWillacv fc Woods! I2OO90 . It is not 
entirely straightforward to relate the column densities 
derived here to those derived from thermo-chemical disk 
models. However, the small scatter in observed molec- 
ular column density ratios, compared to large differ- 
ences between models, suggests that comparisons be- 
tween these data and output of line-generating thermo- 
chemical codes holds great promise for pinpointing the 
most important physical and chemical processes in pro- 
toplanetary disks. 

6. DISCUSSION 

6.1. Completeness of the observed chemistry 

Given the detections of a wide range of abun- 
dant molecules containing three primary elemental con- 
stituents of disks, oxygen, carbon and nitrogen, it is 
possible to compare the observed abundances of these 
species to the total expected elemental abundance, as- 
suming solar values. Obviously, not all molecular species 
have been observed, but what fraction of them are we 
likely to have seen? Using the average column density 
values from Figure Q21 we have accounted for total col- 
umn densities of O ~ 4 x 10 18 cm" 2 , C~lx 10 18 cm -2 , 
and N ~ 2 x 10 15 cm~ 2 , corresponding to ratios of 
O/C ~ 4 and N/O ~ 5 x 10~ 4 . Comparing this to 
the solar ratios of Q /C,„,„„ = 1.6 and N/O = 0.15 

_ ^^^^^^^^^^ bOLG,! ' 

(jFrisch fc Slavinll2003l ). it is apparent that carbon is rel- 
atively underrepresented, while we only detect a very 
small fraction of the available nitrogen. The uncertain- 
ties in the derived abundances, along with the possibility 
that the different species do not probe the same radii, 
probably preclude any strong conclusions to be drawn 
from the derived O/C ratio, except to note that it is con- 
sistent with the observed oxygen-dominated chemistry 
rich in H2O. The missing nitrogen, however, seems to in- 
dicate that a significant chemical reservoir still remains 
to be discovered. Possible carriers of the remaining ni- 
trogen are NH 3 and N 2 . 

NH3 emits strongly near 10/im, so Spitzer-IRS can be 
used to place limits on its column density. With a con- 
servatively low temperature (400 K) and assuming the 
same emitting area as that for the H 2 fits, we find 
upper limits near 10 16 cm~ 2 , or only about 10 x the 
observed column density of HCN. Therefore, NH3 can 
not account for the orders of magnitude of missing N. A 



smaller N deficiency has been noted in comet ary obser - 
vations (IJessbereer fc Kissell Il99lt I Wvckoff et alJH99H) . 
with one explanation being that N 2 is the primary reser- 
voir for N i n the solar nebu la, and is not easily retained 
by comets (jlro et al.1 120 03). Thus, N2 may be a likely 
reservoir for the disks observed in this study as well. 

6.2. The column density of the molecular layer 

The elemental abundances discussed above place strict 
lower limits on the total column density of the warm 
molecular layer at ~1AU. The strongest constraint is 
given by oxygen, even with a significantly unconstrained 
atomic O component. Using the ISM oxygen abundance, 
the thickness of the warm molecular layer must be at 
least N(H) ~7x 10 21 cm" 2 . This is a very g ood match 
to existing models (e.g. iGlassgold et "aDl2009T l . 

6.3. Implications for the carbon deficit in the inner 

solar system 

The terrestrial planets in the solar system have a well- 
known deficit of carbon, by as much as three orders of 
magnitude relative to the elemental abundances in the 
sun. Observations of polluted (by accreting debris) white 
dwarf atmospheres also ind icate simil ar carbon deficits 
in other planetary systems (|Jurall2008h . The high abun- 
dances of carbon observed in the inner disk surfaces ap- 
parently do not allow for a similar carbon deficit; while 
absolute C and O abundances (relative to H) may re- 
main uncertain, a carbon deficit of a factor 10 3 would 
imply that the column densities of warm molecular gas 
are very large, at least AT(H)~ 7V(CO)/(10- 3 a;(C)/) = 
3// x 10 24 cm -2 , where / is the fraction of C accounted 
for in the observed lines, and that a comparatively small 
fraction of the elemental oxygen is accounted for in the 
observations. Both of these requirements seem implau- 
sible in light of recent thermo-chemical models. More 
likely, the observed high abundance of carbon in the gas 
phase indicates that only a small fraction of it will even- 
tually be incorporated into terrestrial planets (however, 
see the discussion of the formation of refractory organics 
in flOJl . 

Recently, iLee et all ((2010T ) proposed that primordial 
carbon grains arc destroyed in the disk surface through 
chemisputtering by free oxygen atoms; essentially they 
burn, producing CO gas. If vertical mixing is efficient, 
these authors suggest that this process can reduce inner 
disk carbon to the observed solar system depletion levels 
as the CO gas dissipates. We do observe the bulk of the 
gas-phase carbon in the form of CO, a significant fraction 
of which could be formed by destruction of carbon grains. 
That there is no obvious observed carbon deficit in the 
gas phase is fully consistent with the picture proposed by 
ILee et afl (|2010l i . If all of the relatively volatile gas-phase 
CO is removed, and assuming there is no carbon left in 
primordial grains, the carbon remaining in CO2, C2H 2 
and HCN is consistent with a carbon deficit of 10 -2 — 
10 -3 . The main uncertainty is likely the abundance of 
CH4, which does not have emission bands in the range of 
the Spitzer IRS high resolution modules. Non-detection 
of free carbon may indicate that it i s not abundant in 
the disk surface (jMathews et al.ll2010f ). 

6.4. Abundant acetylene and aromatic hydrocarbons 
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Many protoplanetary d isks are known to harbor PAHs, 
at lea st in their surfaces ()Habart et al.ll2004l ; iGeers et"aLl 
2007), and some of them have emission spectra entirely 
dominated by the characteris tic mid-infrared PAH emis- 
sion bands (jGeers e t al. 2007). While the detection frac- 
tion is high in HAeBe stars, the relative lack of excit- 
ing UV photons in disks around young solar analogs and 
lower mass stars makes positive detections difficult even 
if PA Hs are actually present in abundance f Gee rs et al.l 
12007( 1 . The observed PAHs are likely analogs to the solar 
nebula precursors of the arom atic component of carbona- 
ceous chondritic meteorites ( jPering fe Ponnamperumal 
Il971f) . which, due to their refractory properties, may rep- 
resent an important route for delivery of organics to the 
surfaces of terrestrial planets. The pathway to PAH pres- 
ence in disks, however, is still debated. 

The formation of aromatics by hydrocarbon combus- 
tion chemistry in the atmosphe res of asymptotic gi- 
ant br anch stars was proposed by iFrenklach fc Feigelsonl 
(1989) and quantified using a detailed kinet ic model. The 
model was applied to the solar nebula by Morgan et al. 
(|1991( 1 to explain the abundance of aromatics in primitive 
chondritic meteorites. Their finding is that the formation 
of PAHs is predicated on (1) the presence of dense gas at 
temperatures of 900-1100K and (2) a high concentration 
of acetylene, C2H2. Their conclusions were, however, 
speculative since it was not known whether the neces- 
sary precursor, C2H2, was indeed present in the early 
solar nebula. 

Our finding that C2H2 is a common constituent of 
protoplanetary disk surfaces at temperatures of 600- 
1200 K therefore has direct implications for the syn- 
thesis of PAHs within typic al protoplaneta r y disk s, in- 
cluding the solar nebula. iMorgan et al.l (|1991( 1 find 
that the PAH yield (fraction of initial carbon atoms 
incorporated into PAH molecules) is a sensitive func- 
tion of pressure as well as the initial C2H2 concentra- 
tion. The most favorable pressures correspond to gas 
densities of -10 14 - 10 15 cm and high abundances of 
C 2 H 2 /H 2 ~10- 6 - 10- 5 , leading to PAH yields (the frac- 
tion of C atoms in the initial gas incorporated into PAHs) 
of order unity. The observed C2H2 abundances are 1-2 
orders of magnitude lower. However, actual PAH yields 
may only be of order a few %, and a significant part of 
the C2H2 may already have been converted to PAHs at 
the typical age of our disk sample (1-3 Myr). Similarly, 
at sufficiently high abundances HCN can polymerize into 
purine and pyrimidine compounds (such as those present 
in nucleic acids), and thus the observation of abundant 
C2H2 and HCN at 1000 K supports the possibility of a 
complex organic/ "soot ring" chemistry in protoplanetary 
disks, which may be revealed by high resolution imaging, 
including mid-IR interferometry. 



Alternatively, it has been shown that starting with 
abundant PAH's, chemical reactions with H, OH and O 
tend to destroy PAHs in the terre strial planet-forming 
regions of disks (|Kress et all 12010( 1 . One by-product of 
this process is C2H2, with a maximum C2H2/CO ratio 
of ~ 10~ 2 . In either case, the region of the disk sur- 
face traced by the mid-infrared molecular line forest may 
indeed be an important site for understanding the syn- 
thesis and/or destruction of the carbonaceous, aromatic 
component of chondrites. 

6.5. Conclusion 

This work presents an overview of the rich near- 
and mid-IR molecular emission observed from proto- 
planetary disks. What physical mechanisms account 
for the presence or absence of molecular lines, and for 
differences in molecular line ratios from source to source? 
What physical process are most crucial for determining 
chemical abundances in the inner disk, and is there in- 
heritance of information from the cloud it formed out of? 
We have begun to answer these questions, showing here 
and in Paper I the importance of photochemistry, disk 
heating, grain settling and various chemical processes in 
the production of molecular line emission. However, it 
is also clear that there are variations within the sample 
that have yet to be explained. In the near future, these 
questions will continue to be explored, with radiative 
transfer codes and thermo-chemical disk models. We 
also expect that the data presented here will provide a 
rich database for use with observations obtained both in 
the near future and beyond, with facilities such as the 
Herschel Space Observatory, ALMA, SOFIA, JWST, 
and the emerging class of 'extremely large' telescopes. 
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